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Abstract 

We consider the approach to self-similarity (or dynamical scal- 
ing) in Smoluchowski's coagulation equations for the solvable kernels 
K(x, y) = 2,x -\- y and xy. We prove the uniform convergence of den- 
sities to the self-similar solution with exponential tails under the reg- 
ularity hypothesis that a suitable moment have an integrable Fourier 
transform. For the discrete equations we prove uniform convergence 
under optimal moment hypotheses. Our results are completely anal- 
ogous to classical local convergence theorems for the normal law in 
probability theory. The proofs rely on the Fourier inversion formula 
and the solution by the method of characteristics for the Laplace trans- 
form. 

1 Introduction 

Smoluchowski's coagulation equation 

dtn{t,x) = - K{x-y,y)n{t,x-y)n{t,y)dy- K{x,y)n{t, x)n{t,y)dy, 
^ Jo Jo 

(1.1) 

is a widely studied mean-field model for cluster growth [4, 8, 17]. We study 
the evolution of n{t,x), the number of clusters of mass x per unit volume 
at time t, which coalesce by binary collisions with a symmetric rate kernel 
K(x,y). Equation (1.1) has been used as a model of cluster growth in a 

'^Department of Mathematics, University of Wisconsin, Madison WI 53706. Current ad- 
dress: Division of Applied Mathematics, Brown University, Providence, RI 02912. Email: 
menon@dam.brown.edu 

^Department of Mathematics & Institute for Physical Science and Technology, Univer- 
sity of Maryland, College Park MD 20742. Email: rlp@math.umd.edu 



1 



surprisingly diverse range of fields such as physical chemistry, astrophysics, 
and population dynamics (see [4] for a review of applications). In addition, 
over the past few years a rich mathematical theory has been developed for 
these equations. Aldous [1] provides an excellent introduction. 

Many kernels in applications are homogeneous, that is K{ax,ay) = 
a'^K{x,y), x,y,a > 0, for some exponent 7 [4]. A mathematical problem of 
scientific interest is to study self-similar or dynamical scaling behavior for 
homogeneous kernels. There are no general mathematical results for this 
problem despite an extensive scientific literature (especially formal asymp- 
totics and numerics [12, 13, 18]). It is known that 7 plays a crucial role. 
On physical grounds, we expect solutions to (1.1) to conserve the total mass 

xn{t, x)dx. When K(x,y) < 1 + x + y (corresponding to < 7 < 1), 
mass-conserving solutions exist globally in time under suitable moment hy- 
potheses on initial data [5]. It is then typical in applications to assert that 
the solutions approach "scaling form" [13, 18], but there is no rigorous math- 
ematical justification for this in general. 

For a large class of kernels satisfying (xy)"'/'^ < K{x, y) with 1 < 7 < 2, 
it is known that there is no solution that preserves mass for all time. This 
breakdown phenomenon is known as gelation. It was first demonstrated by 
McLeod [14] with an explicit solution for the kernel K = xy. A general 
result using only the growth of the kernel was proved probabilistically by 
Jeon [9] (see also [6] for a simple analytical proof). It is natTiral to ask 
whether the blow-up is self-similar, but there are no general results on this 
problem yet. 

There are a number of results, however, for the 'solvable' kernels K = 2, 
X + y and xy (sec [15] and references therein; also see [13]). A remarkable 
feature of these kernels is that the problem of dynamical scaling can be un- 
derstood quite deeply by analogy with classical limit theorems in probability 
theory. For example, an analog to the classical Levy-Khintchine representa- 
tion for infinitely divisible laws was proved by Bertoin [2] for eternal solutions 
to Smoluchowski's equation with kernel K = x + y. Eternal solutions are 
defined for all t G (—00,00), meaning they model coagulation processes 'in- 
finitely divisible' under Smoluchowski dynamics. Later, we proved [15] that 
the domains of attraction of self-similar solutions (in the sense of weak con- 
vergence of measures) can be characterized by almost power-law behavior 
of the tails of the initial size distribution. This is analogous to the char- 
acterization of the weak domains of attraction of the Levy stable laws [7]. 
An essential component in both proofs is a simple solution formula for the 
Laplace transform of n that is widely known [4]. These results may be used 
as a basis for refined convergence theorems, as we now explain. 
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A general theme in probabilistic limit theorems is the interplay between 
moment and regularity hypotheses and the topology of convergence. In this 
article, we develop one aspect of this idea. Under stronger regularity hy- 
potheses, the weak convergence results of [15] will be strengthened to obtain 
uniform convergence of densities using the Fourier transform. This method 
is classical in probability theory and is used to prove uniform convergence 
of densities in the central limit theorem [7, XV. 5. 2]. Feller's argument in 
[7] is simple and robust, and our main contribution is to show that it ex- 
tends naturally to Smoluchowski's equation. The key new idea is to use the 
method of characteristics in the right half of the complex plane to obtain 
strong decay estimates on the Laplace transform. A broader contribution 
of this work and [15] is to show that the analytical methods used to prove 
classical limit theorems in probability apply to a wider range of problems 
involving scaling phenomenon for integral equations of convolution type. 

Let us briefly connect our results to previous work: the only uniform 
convergence theorems in the literature are that of Kreer and Penrose for the 
kernel K = 2 [11], and closely connected work of daCosta [3]. In this article, 
for K = 2 and a; -|- y we present theorems on uniform convergence to the 
self-similar solutions with exponential tails for the continuous and discrete 
Smoluchowski equations. For K = xy, we prove uniform convergence of 
densities to self-similar form as t approaches the gelation time T^^i- For 
K = 2, wc strengthen the result of Kreer and Penrose and simplify the proof. 
Their decay hypothesis on the initial data (no(x) < Ce^""^) is weakened to 
an (almost) optimal moment hypothesis, and their regularity hypothesis 
(no G C^) is weakened to a little bit more than continuity. For K = x + y 
the convergence theorem is new. Study of the kernel K = xy is reduced to 
K = X + y hy a well-known change of variables [4] . Uniform convergence to 
the self-similar solutions with "fat" or "heavy" tails is a more delicate issue, 
and will not be considered here. 

Our uniform convergence theorems may be stated in a unified manner as 
follows for the continuous Smoluchowski equations with kernels K{x, y) = 2, 
X + y and xy, corresponding to 7 = 0, 1, 2 respectively. Presuming the 7-th 
and (7 -I- l)-st moments are finite, we may scale x and n so both moments 
are initially 1. For the multiplicative kernel this ensures that the gelation 
time Tgci = 1. Let = 00 for 7 = 0, 1, = Tgd = 1 for 7 = 2. The 
self-similar solutions with exponential tails are explicitly given by [15] 




(1.2) 
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where for x >0, 

n*fi{x) = e~^, xn*,i(x) = x^n*,2(^) = -^x~-^/^e~^/^, (1.3) 

V 27r 

and 

mo(t) = r\ mi(t) = l, m2(t) = (l-t)-\ (1.4) 

Ao(0 = i, Ai(0 = e2*, X2{t) = il-t)-\ (1.5) 

Our sufficient conditions for uniform convergence to these self-similar so- 
lutions for the continuous Smoluchowski equations are summarized by the 
following result. 

Theorem 1.1. Let no > 0, x'^nQ{x)dx = x^^'^nQ(x)dx = 1. Assume 
that the Fourier transform of x^^^hq is integrable, and let n{t, x) be the 
solution to Smoluchowski' s equation with initial data no(x) and K = 2, 
X + y or xy, for 7 = 0, 1 or 2. Then the rescaled solution 

h{t,x) = ^ n(t, xXj{t)) 

satisfies 

lim sup x^'^'^ |n(t, x) — n*^^(x)| = 0. 

It has been traditional to treat the discrete Smoluchowski equations sepa- 
rately from the continuous equations. Yet, within the framework of measure 
valued solutions [15, 16], the discrete Smoluchowski equations simply corre- 
spond to the special case of a lattice distribution, a measure valued solution 
supported on the lattice hN and taking the form = ^^i(i)^fei(-^;), 
where 6fii{x) is a Dirac delta at hi. If /i is maximal we call a lattice 
measure with span h. The coefficients ni satisfy the discrete Smoluchowski 
equations 



^ i—i 00 
9tni{t) = 2 '^'-JJ^i -J - m '^l,3Mt)nj{t), (1.6) 



where nij = K(lh,jh). Physically, this case is of importance, since some 
mass aggregation processes (e.g., polymerization) have a fundamental unit 
of mass (e.g., a monomer). The uniform convergence theorems for the con- 
tinuous Smoluchowski equations have a natural extension to this case. 
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Theorem 1.2. Let > be a lattice measure with span h such that 
Jq°° x't'i^o (dx) = x^~^'^i'o{dx) = 1. Then with 



A-y(i)' h m-y{t) 



we have 



lim sup 

[1+7 



ni{t) - n^,^{i) 



0. 



Let us comment on the hypotheses in Theorems 1.1 and 1.2. The moment 

hypotheses in both theorems are essentially the same. x^uo{dx) = 1 is 
the natural hypothesis for existence and uniqueness of solutions [15]. The 
other moment condition x^~^^i'o(dx) = 1 is of a different character. It 
implies that no or is in the weak domain of attraction of the self-similar 
solution with exponential tail, under a rescaling n{t, x) — >■ h{t, x) that fixes 
both moments 



Jo 



1 /•OO 

x'^n{i, x)dx = / x'^'^^h{t, x)dx = 1 for all £ > 0. 
Jo 



The hypothesis that the (7 + l)-st moment is finite is almost optimal. The 
weak domain of attraction under a broader class of rescalings is a bit bigger, 
as it allows for a weak divergence x^'^^v'Q{dx) ~ L{y) as y ^ 00 for a 
slowly varying function L(y) [15]. Thus, Theorem 1.2 shows that within the 
class of lattice measures, the weak convergence of measures almost implies 
uniform convergence of the coefficients. 

Theorem 1.1 requires an additional hypothesis on integrability of a suit- 
able Fourier transform. This is a regularity hypothesis that is the analog of 
the hypothesis for uniform convergence to the normal law used by Feller [7]. 
One may heuristically understand the role of regularity as follows. Equation 
(1.1) is hyperbolic and discontinuities in the initial data persist for all finite 
times. On the other hand, the self-similar solutions in (1.3) are analytic. 
Thus, one expects some regularity on the initial data is necessary to obtain 
uniform convergence to a self-similar solution. Loosely speaking, regularity 
of the initial data no{x) translates into a decay hypothesis on its Fourier 
transform. We need only the weak decay implied by integrability. 

We do not know if this assumption is optimal, or if it may be weakened 
further. We briefly comment on this issue here; it will not be considered 
in the rest of the paper. The space of functions with integrable Fourier 
transforms is of great interest in harmonic analysis. Precisely, for / € 
let F be its Fourier transform. Then the space 

^(R) = {/ G L^(M)|F G 
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is a closed subalgebra of L-'^(M) known as the Wiener algebra [10]. Inte- 
grability of F implies that / is continuous. But it also implies more. It 
is known that functions in A{M.) possess some delicate regularity proper- 
ties. For example, a function in A(R) has a logarithmic modulus of con- 
tinuity in a neighborhood where it is monotonic. It is definitely not ob- 
vious whether this regularity is truly necessary to obtain uniform conver- 
gence. If VQ{ik) = e^'^''^x^^'^nQ{x)dx is integrable it also follows that 
vq G i/^(M) n ^(M), since vq is the boundary limit of an analytic function 
(the Laplace transform of x^'^'^uq). Here denotes the classical Hardy 
space. This is turn means that vq has some hidden regularity and integra- 
bility properties. It is worth remarking that the precise characterization of 
A(R) remains an outstanding open problem in harmonic analysis (though 
several sufficient conditions are known, see [10]). 

2 Uniform convergence of densities for the con- 
stant kernel K = 2 

2.1 Evolution of the Laplace transform 

Let C+ = {z G C I Re z > 0} and C+ = {z G C | Re z > 0}. We let 



denote the Laplace transform of the number density n. We take the Laplace 
transform of (1.1) with K = 2, and its limit as z — to see that N{t,z) 
solves 



Without loss of generality, we may suppose that the initial time t = 1. We 
will always assume that the initial data is normalized such that 



If the initial number of clusters, n(l, x)dx, and the mass, xn{l, x)dx, 
are finite, we may always assume (2.2) holds after rescaling x and n. We 
solve the second equation in (2.1) to see that the total number of clusters 
decreases according to 




dtN = N^- 2N{t, 0)N, dtN{t, 0) = -N{t, Of. 



(2.1) 




(2.2) 




-1 



t > 1. 



(2.3) 
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We hold z fixed and integrate (2.1) in t to obtain the solution 

Mt.) = l (2 4) 

tt(l-iV(l,z)) + iV(l,z)- ^^-^^ 

The evolution preserves mass. Indeed, if we differentiate (2.4) with respect 
to z, we find 

/•cx) /•oo 

/ xn{t, X) dx = -d^N{t, 0) = -d^N{l, 0) = / xn{l, x) dx = 1. (2.5) 
Jo Jo 

2.2 Approach to self-similarity 

A special case of the weak convergence result of [15], also given by Leyvraz [13], 
is obtained as follows: Observe that for each fixed s G C_)_ equations (2.3), 
(2.4), and (2.5) imply 

tmt St-') - .J- (2 6) 

It is classical that the pointwisc convergence of Laplace transforms is equiv- 
alent to weak convergence of measures [7, XIII. 1.2a]. Thus, (2.6) implies 
that rescaled solutions to Smoluchowski's equations converge weakly. Let 
us be more precise about the rescaling. We define the similarity variables 

x 

T = logt, x = — = e~'^x, s = tz = e^z, (2.7) 
and the rescaled number distribution, 

n(r, x) = e^^n(e^, e'x) = t^n{t, x). (2.8) 
Observe that this rescaling preserves both total number and mass, that is 

/•oo /•oo 

/ n{T,x)dx= / xn{T,x)dx = l, r > 0. (2-9) 
Jo Jo 

We denote the Laplace transform of n(r, x) by 

/•oo 

i(r, s)= e-'^n{T, x) dx = e^Ar(e^, se'^) = tN(t, z). (2.10) 
Jo 

In these variables, the pointwise convergence of (2.6) takes the simple form 

lim u(t,s) = -— — =:u^o{s), s G C+, (2-11) 

T^OO 1 + s 
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where u*fi{s) denotes the Laplace transform of 

n*,Qix) = e"*, x>0, (2.12) 
the profile for the self-similar solution in (1.2). Now, (2.11) is equivalent to 

n(T, x) dx — > n^fl{x) dx. 

as r ^ oo, in the sense of weak convergence of measures. 

Our goal is to strengthen this to uniform convergence in both continuous 
and discrete cases, under appropriate hypotheses on initial data. For the 
continuous Smoluwchowski equation (1.1) we prove 

Theorem 2.1. Let n(l,a;) > 0, n{l,x)dx = xn{l,x) dx = 1. As- 
sume that the Fourier transform of xn{l,x) is integrable. Then in terms of 
the rescaling in (2.7)-(2.8) we have 

lim sup x|n(r, x) — n* o(^)| = Oj (2-13) 

where r7,*^o(^) = is the similarity profile in (2.12). 

The proof of this theorem extends to treat uniform convergence of coeffi- 
cients for solutions of the discrete equations (1.6) under only the hypothesis 
that the zeroth and first moments are finite; see Theorem 2.2 below. 

Observe that we prove uniform convergence of the weighted densities 
xn(r, x). The reason can be ascribed to use of the Fourier-Laplace inversion 
formula. We cannot apply the inversion formula directly to u*^o a-s it is not 
integrable on the imaginary axis (|'U*,o(ifc)| ~ 1^1^^ \k\ —>■ oo). The slow 
decay of the Fourier transform is caused by the jump discontinuity at ar = 0, 
since h^^o{x) = for x < 0. In order to gain a uniform convergence result, 
we smooth this discontinuity and consider the mass density xn. Its Laplace 
transform we denote by 

/•oo 

v{t, s) = -dsu{T, s) = / e~*^fn(r,x)dx. (2-14) 
Jo 

Differentiating (2.11), we obtain a corresponding self-similar profile, with 
^*,o(s) := T^^^- \v*,o{ik)\ = kGR. (2.15) 
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2.3 Evolution on characteristics 



The explicit solution for u{t, s) and v{t, s) can be obtained directly by sub- 
stituting (2.10) into (2.4). But we rederive the solution to make explicit 
the geometric idea underlying the proof of Theorem 2.1. The same ideas 
underlie the proof of Theorem 3.1 for the additive kernel and arc more easily 
understood here. We use the change of variables (2.7) and (2.10) in (2.1), 
and the conservation of moments in (2.9), to obtain the equation of evolution 
for u: 

drU + sdgU = —u{\ — u). (2-16) 

The solution of equation (2.16) may be described by the method of charac- 
teristics. A characteristic curve s{t;tq,sq) is the solution to 

ds 

— = s, s(r;ro,so) = So e C+. (2.17) 

Explicitly, 

s(r;ro,so) = e^-^°so- (2.18) 

Equation (2.17) is an autonomous differential equation in C+, and may be 
thought of geometrically. For fixed sq G C+ the trajectory of the character- 
istic curve s(r;ro,so),r G M, is a ray in C+ emanating from the origin. In 
particular, the imaginary axis is invariant under the flow of (2.17). Equation 
(2.18) shows that the characteristics expand uniformly outward at the rate 
. Along characteristics we have 

^ = -u{l-u), (2.19) 
which may be integrated to obtain the solution 

We need to estimate the decay of the derivative v = —dgU. Differentiating 
equation (2.16), we see that on characteristics the derivative solves 

^ = -2{l-u)v. (2.21) 
ar 

We integrate (2.21) using (2.20) to find 

vir,s) = ^-o^^^y-''-'' ^ . (2.22) 

(l-n(To,so)(l-e-(--o)))2 
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For r > To we may take absolute values in (2.20) and (2.22) to obtain the 
decay estimates 

, , |'u(To,So)|e-(''-''o) 

\uiT,s)\ < / -, r-, (2.23) 

' ^ ' ^' - l-|n(ro,so)|(l-e-(^-^o))' ^ ^ 

and 

Hr,s)\< Kro,.o)|e--(-°)_ 

' ^ ' {l-\u{ro,so)\{l-e-(r-ro))f - il-\uiTo,so)\f 

2.4 Proof of Theorem 2.1 

1. We use the Fourier-Laplace inversion formula 

x{h{T,x) - ra*,o(^)) = I e"^"" (v{T,ik) - v^^o{ik)) dk. (2.25) 
Thus, in order to prove (2.13) it suffices to show 

lim / \v{T,ik)-v*o{ik)\dk = 0. (2.26) 

2. Let e G (0, ^) and put R = e~^. We will prove (2.26) by estimating the 
integral separately in three regions: |A;| < -R, -R < |A;| < Re^~^ and Re^"^ < 
\k\, for T >T where T > will be chosen sufficiently large, depending on e 
and the initial data vq. This is essentially the same decomposition used in 
the proof of uniform convergence in the central limit theorem by Feller [7, 
XV.5.2]. The main new idea here is the use of the decay estimates (2.24) 
and the method of characteristics in the regions where R < \k\. 

3. \k\ < R: Recall that the pointwise convergence of Laplace trans- 
forms (2.11) is equivalent to h{T,x)dx h.tfi{x)dx in the sense of weak 
convergence of measures. Combined with (2.9) this also implies that the 
mass measures xh{T, x) dx converge weakly to xh^^Q{x) dx as r ^ oo. But 
this implies v{T,ik) converges to v*,o(^fc) uniformly for \k\ < R [7, XV.3.2]. 
Therefore, 

lim / \v{T,ik) -v^o{ik)\dk = 0. (2.27) 

4- It remains to consider \k\ > R. It is sufficient to consider only k > R, 
since \v{T,ik)\ = \v{t, —ik)\. We will control v{T,ik) and v^^ separately: 

/ \v{T,ik) — v^fi{ik)\ dk < / \v{T,ik)\dk+ / \v^^o{ik)\dk. 

JR ' JR JR ' 
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But \v*fi{ik)\ = (1 + by (2.15), so that 

\v*fl{ik)\ dk < = e. 



PCX 

JR 



In the rest of the proof we estimate \v{t, ik)\dk. 

5. Since u{T,ik) u.tfi{ik) and v{T,ik) v^^q as r — > oo for each real 
k, using (2.11) and (2.15) we may choose T > such that 

sup \u{t, iR)\ < = e, sup \v{t, iR)\ < R^^. (2.28) 

T>T T>T 

6. R < k < Re^"^: The control obtained from (2.28) propagates out- 
wards along characteristics as r increases. Precisely, whenever t > T, for 
any k such that R < k < Re^~'^ we have ik = e^~'^'^iR where tq > T. 
By (2.18) this means that ik = s(r;ro,so) with so = iR- Then the decay 
estimate (2.24) and the boundary control (2.28) imply 



|z;(To,ii?)|e-2(--o) 1 (R 



2 



" - {l-\u{To,iR)\f - (1-e)' \kj - 
Integrating this estimate we obtain 

i-Re^~^ roo 

/ \v{t, ik)\ dk< 4fe-2 dk = AR'^ = Ae. 

JR JR 

7. Re'-'^ < k: For brevity, let R = Re''^ . With uo{s) := u{0,s), 
vo{s) := v{0, s), we use (2.24) and (2.18) with tq = to obtain 



" Ar,ik)\dk < e-^^ r ^ dk 

Re-^ JRe^ (1 - \uo{ike ^)|) 



/ odk < sup Ty e ko Li, 

Jr {l-\uo{ik>)\f - l|fc,|>R(l-|no(ifcO|)V 



W\>R 

where k' = ke~'^ . Since \uQ[ik')\ < 1 for k' ^ and UQ{ik') ^ as ^ oo 
by the Riemann-Lebesgue lemma, we have sup|^,|>^(l — |uo(ifc')|)~^ < oo. 

8. Putting together the estimates we have obtained, it follows that for r 
sufHciently large, the integral in (2.26) is less than 12£. This completes the 
proof. 
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2.5 The discrete Smoluchowski equations 



We consider measure solutions of the form ut = Yl'^i''^l{'t)^hl{x), where 
6hi{x) denotes a Dirac mass at hi. To avoid redundancy, we always assume 
that h is the span of the lattice, that is, the maximal h > so that all 
initial clusters, and thus clusters at any time t > 0, are concentrated on 
hN. We will call ut a lattice measure with span h. Notice that if the initial 
number of clusters and the mass are finite, by rescaling n; and h we may 
assume that iyi{dx) = xui{dx) = 1. Under these conditions, the 
weak convergence theorem of [15] asserts that \\mt^^tN{t^s/t) = m*,o('S)- 
We show that this theorem may be strengthened by use of Fourier series. 
The Fourier transform of vt is the Fourier series 

N{t,ik) = ^n^(^)e-^'''^ G M, 
which has minimal period 2Tx/h. T\msni{t) = {h/ 27r) J^l^^^e'^''''N{t,ik) dk, 



or 



t ni{t) = 7^ / ex.'p{ilhke~'^)u{T,ik) dk, 

in similarity variables from (2.10). We integrate by parts and let 

i=lhe~'' = lht~^, hi{t) = h'^t^ni{t) 

to obtain 

e''^v{T, ik) dk. 



1 r-Ke^ jh 

ihi{t) = tlni{t) = — / 



(2.30) 



(2.31) 



(2.32) 



-Tve'^ /h 

As in Theorem 2.1 we expect the right hand side to converge to ln^fi{l) as 
r — ^ oo, indeed uniformly for I G ht~^N. 

Theorem 2.2. Let > be a lattice measure with span h such that 
Jo° ^i{dx) = xvi{dx) = 1. Then with the scaling (2.31) we have 



lim sup I hi{t) — n^fi{l) 



t— »oo 



(2.33) 



Proof. By (2.32) and the continuous Fourier inversion formulas it suffices to 
show that 



lim sup 

/>o 



t 



-ilk. 



-J 

Jk 



v{T,ik)dk- I e^^^v^fi{ik)dk 



= 0. 
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As earlier it suffices to consider A; > 0. Tlie integrals 



/n pne' poo 

\v{T,ik) — V:^^o{ik)\dk, / \v{T,ik)\dk, / \v^fl{ik)\dk, 
-R ' Jr Jr ' 



f 

Jr, 



with R = Re~^ , are controlled exactly as in the proof of Theorem 2.1. It 
only remains to estimate the integral of |u(r, ik)\ over the region Re^ < k < 
Tve^ /h. We assume that 7r/h > R, for otherwise there is nothing to prove. 
But then by the formula (2.18), the uniform decay estimate (2.24), and the 
change of variables k' = fee^^, we have 

i-Tre^/h i-n/h \v^^(^k')\ 

\v{T,ik)\dk<e-^ ^^^^^ .dfc'. 

Since the domain of integration is finite, it suffices to show that the integrand 
is uniformly bounded in time. Since |uo(ifc)| < 1, it is only necessary to 
control the denominator. But uo{ik) = X^/gj^ "KO)^"'"^'^ with n/(0) > 0. 
Therefore, \uo{ik)\ < 1, and [7, XV.1.4] yields that 

uo{ik) = 1 if and only if k = — ; — , m G Z. 

h 

In particular, we have the strict inequality 



min 11 — uo{ik)\ > S > 0. 

k&[R,l] 



Therefore, 



|l - UQ{ik){l - e-^)\ > |1 - UQ{ik)\ - \uQ{ik)\e-^ > 5 - e'^ > - 
for sufficiently large r. Thus, 

/ \v{T,ik)\dk < —e-\ 



□ 



3 Uniform convergence of densities for the addi- 
tive kernel 

3.1 Rescaling and approach to self-similarity 

In this section we prove the analogs of Theorems 2.1 and 2.2 for the additive 
kernel. The essential geometric ideas of the proof are similar to the previous 
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section. However, the trajectories of the characteristic curves s{t;to,so) in 
the complex plane are no longer rays, and the proofs require more careful 
analysis. As earlier, we will work with the explicit solution formula for an 
appropriate Laplace transform. For z G C+ we define 

/•oo 

z)= (1 - e-^"") n{t, x)dx. (3.1) 

We observe that 1 — e~^^ = zx + 0{z'^x'^) as x ^ 0. We use $ instead of the 
standard Laplace transform of n because the latter may not be well-defined: 
E.g., the similarity profile n*,i in (1.3) satisfies n*,i(x) ~ Cx~^l'^ as x — > 0. 
More generally, one needs the initial data to have only a finite first moment 
for existence and uniqueness of a solution to (1.1) in the case of the additive 
kernel [15]. A deeper reason for this choice of variables (and notation) is 
probabilistic: (3.1) is the Levy-Khintchine formula for the Laplace exponent 
of a subordinator with no drift [2]. We will always assume that the initial 
data no satisfies the moment conditions 

/•oo /"OO 

/ xnQ{x)dx = 1, / x^nQ{x)dx = 1. (3-2) 
Jo Jo 

We substitute (3.1) in (1.1) and use (3.2) to see that ^{t,z) solves the 
equation 

POO 

dt^ - <^d,^ = $(0, z)= (1 - e-''')no{x)dx. (3.3) 

Jo 

As shown in [15] by the method of characteristics, (3.3) has a unique 
solution for z > 0, t > which is analytic with derivative completely 
monotone in z and satisfying dz^{t,0) = 1 for all t. For each t > then, 
dz^{t, •) is the Laplace transform of a probability measure, so its domain 
contains C+ and (3.3) holds by analytic continuation for z G C+, t > 0. 

In contrast with (2.4) it is not obvious that a suitable rescaling will lead 
to convergence to self-similar form. This point is discussed in [15, Sec. 7], and 
we refer the reader to that article for motivation for the following change of 
variables. We define the similarity variables 

X = xe""^*, s = ze^\ (3.4) 

and the rescaled number density 

n(i, x) = e'^^n{t, xe^*) = e^*n(t, x). (3.5) 
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We also define the rescaled Laplace transforms 

roc 

ip{t, s) = e2*$(i, e-2*s) = / (1 - e-'^)n{t, x)dx. (3.6) 

Jo 

Part of the motivation for the rescahng (3.4) and (3.5) is that this choice 
preserves both moment conditions in (3.2). That is, we have 

/•oo /"OO 

/ xn{t,x)dx = I x'^n{t,x)dx = 1, t>0. (3.7) 
Jo Jo 

This should be compared with (2.9) for the constant kernel. The mass 
measure plays the same role here as the number measure did for K = 2. 
Thus, we denote its Laplace transform by the same letter, and let 

/•oo 

u{t, s) = ds(p{t, s) = / e~^'-^xn{t,x)dx. (3.8) 
Jo 

By Theorem 7.1 in [15] (also see [13, Appendix G]), the assumptions in (3.2) 
imply that the rescaled mass measures converge to the similarity profile, with 

xfi{t,x)dx ^ xh*^i{x)dx = —=x~^^'^e~^^'^dx, t — > oo (3.9) 

V 27r 

in the sense of weak convergence of measures. It then follows from [7, 
XIII.1.2] that (3.9) is equivalent to 

lim u(t, s) = —j^=^= =: i(s). s G C+, (3.10) 
*-^oo Vl + 2s 

Our goal is to strengthen (3.9) to uniform convergence of densities for 
(1.1) and uniform convergence of coefficients for (1.6). For the continuous 
Smoluchowski equations we prove 

Theorem 3.1. Suppose no{x) > 0, xno{x)dx = x^no{x)dx = 1. 
Suppose also that the Fourier transform ofx^riQ is integrable. Then in terms 
of the resettling (3.4)-(3.5) we have 

lim sup x^|n(t, x) — i(x)| = 0, (3-11) 

x>0 

where n*,i(x) is the similarity profile defined in (1-3). 
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Once Theorem 3.1 is established, it is relatively straightforward to ob- 
tain the analogous result for the discrete Smoluchowski equations; see The- 
orem 3.6 below. Thus, most of our effort is devoted to Theorem 3.1. 

Observe that we prove uniform convergence of the weighted density 
x^h{t,x). As in the previous section, this is because Theorem 3.1 is proved 
using the Fourier-Laplace inversion formula. Since |MH,,i(ifc)| ~ as 
\k\ oo, tt^^i is not integrable on the imaginary axis. This divergence is 
due to the fact that n*,i(a;) = for a; < and xn*^i(a;) ~ Cx~^/^ as x ^ 0+. 
As earlier, we resolve the situation by considering the transform of the next 
moment. Let 

/•oo 

v{t,s) = -dsu{t,s) = / e~^^x^n{t,x)dx, s G C+. (3.12) 
Jo 

We integrate and differentiate (3.10) to obtain 

<p^,i{s) = - 1, v^,i{s) = {l + 2s)-^/^, seC+. (3.13) 

3.2 Characteristics and estimates 

The equations of evolution for (p and u are 

dt(p + {2s - <p)dsip = ip, (3.14) 
dtu + {2s - (p)dsu = -u{l-u). (3.15) 

In what follows, we first derive solution formulas to (3.14) by the method 
of characteristics. We then show that the solution map for the characteristic 
equation is never degenerate and that characteristics flow out of the right 
half into the left half of the complex plane as t increases. For most parts 
of our analysis, it will suffice to study characteristics in the right half plane 
only. But for one part, we need to study characteristics that start in the 
right half plane but move into the left half plane. 

We use the notation s{t;to, so) to denote the solution to 

ds 

— = 2s-ip, s{to;to,so) = sq. (3.16) 
Along the characteristic curve s{t;to,so) we have 

1 = ,. and | = -m-.). (3.17) 
We integrate (3.17) to obtain 

^{t, s) = e*-*V(to, so), u{t, s) = — -. (3.18) 

1 - u{to,so){l -e *oi) 
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We now substitute for ip{t,s) from (3.18) in (3.16) and integrate to obtain 
the explicit solution 

g-2(t-to)^(i. to, so) = SO - ^{to, so){l - e-(*-*°)). (3.19) 

This equation can also be rewritten in two other useful forms, namely 

e-2{t-to) _ = - ^{to, so)) , (3.20) 

and 

V>{t,s) ^ (y(to,go)/so)e-(*-*°) 

l-(<^(to,5o)Ao)(l-e-(*-*o))" ^ ■ ' 

The method of characteristics also yields an explicit solution for v{t, s). We 
differentiate equation (3.15) to obtain 

^ = -3(1 - u)v. (3.22) 
We substitute for u from (3.18) and integrate (3.22) to obtain 

"('o>'-°)<=-'"-'°' J. (3.23) 

Let </9o('S) := ¥'(0, s) and similarly mo(s) := u{0,s), vo{s) := u(0, s). 
Since u = dgip and ip{t, 0) = 0, the moment conditions (3.2) and the identity 
ipo{s)/s = uo{ts) dr imply 

\uo{s)\ < 1, \vo{s)\ < 1, \Ms)\ < \s\, s e C+. (3.24) 

These inequalities are strict for s 7^ because xno{x)dx is not a lattice 
measure [7, XV. 1.4]. Taking to = at first, for t > to we take absolute 
values in (3.18) and (3.23) to see that |u| and l^l decay along characteristics 
according to 

I ^| ^ |^(to,so)|e-(*-*») 

1 - \u{to,so)\[l -e *oj) 

\v{to,so)\e-^^'-'° ^ 
{l-\u{to,so)\)^ 

From (3.25) and the fact that |tto(so)| < 1 for so / 0, and a similar 
estimate using (3.21) and \ipo{so)/so\ < 1, it follows 



Kt,s)i< 3 • (3.26) 



\u{t,s)\<l, \(p{t,s)/s\ <1, t>0, s^O. (3.27) 
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Then (3.25) and (3.26) hold also for any to > 0, if t > to- 

Let us also note the uniform outward growth of characteristics implied 
by (3.27). Using (3.27) together with (3.16) we obtain 



N<M<3|4 



(3.28) 



Thus, |so|e^*^*"^ < |s| < e'^''*~*"^|so|- We will refine this crude estimate in 
the proof of Theorem 3.1, but we note here that |s(t;io!So)| is a strictly 
increasing function of t. 

In addition to the decay along characteristics, we will need the following 
uniform Riemann-Lebesgue lemma. Let Cji = {s E \ \s\ = R} denote 
the semicircle of radius R in the right half plane. 



Lemma 3.2. Let g{x) e L^{0,oo) and G{s) = e ^^g{x)dx. Then 



lim sup \G{s)\ = 0. 



(3.29) 



Proof. Let e > 0. We choose a step function = Yl^=i'^k'^[ak,bk] that 
Wq ~ ffellLi < But then, ||e~*^(5' — ^e)!!^! < s. Therefore, for s G C+, 



\G{s)\<e + 



''ge{x)dx 



e + 



^ rbk 
fe=i -^"fc 



'dx 



\s\ 



□ 



We apply this lemma and (3.7) to g{x) = x^h{t,x) for j = 1,2 to infer 
that for every t > 0, as |s| — oo with Res > we have 



\u{t,s)\^0, \v{t,s)\^0, 



^{t, 8) 



0. 



(3.30) 



3.3 Geometry of the chciracteristic map in the complex plane 

In this subsection, we study the solution formula (3.19). Our goal is to 
delineate some key properties of the map sq i— > s{t; tQ, sq) for t, to > 0. 

Let C_|_ denote the open right half plane. We let fij denote the image 
of C+ under the map sq h- > s{t;0,SQ), and let Tf denote the image of the 
imaginary axis under the same map. We aim to prove the following. 

Lemma 3.3. (i) For any t > Q, Tt is a curve that passes through the 
origin but otherwise lies in the open left half plane. On Tf, Res is a 
function o/Ims. 
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(a) fit is the component of the complex plane to the right ofFf. Conse- 
quently Tf = d^t o-'f'd fit D C+ \ {0}. 

(Hi) Whenever ti >tQ>Q, the map sq ^ s\ = s(ti;fo)So) is one to one 
from fito onto fi^j. It is on fifj, and analytic in fi^y. The inverse 
map is given by si i— > sq = s{to;ti, si), and is on fi^^ and analytic 
in ■ 

(iv) Whenever ti > and si G C+, the backward characteristic curve 
s(to; ii) -si), to G [0,ti], lies mC+. 

Proof. We first establish part (iii), taking to = at first. Since x^uq is 
integrable, vq{s) is continuous in C+ and analytic for Res > 0. It follows 
by a standard dominated convergence argument that uq is and (po is 
in C+, and these functions are analytic in C+. Prom (3.19) we see that the 
map So 1-^ s(t;0,so) is analytic in C-|- and on C^. (meaning derivatives 
up to second order extend continuously to C+). 

We next claim that this map is one to one. The proof relies on the fact 
that (^0 is contractive, with 

Wa{so) - 'Po{so)\<\sq- so\, so,so,gC+. (3.31) 

This holds because |5sV3o(s)| < 1 for s G C+ as an immediate consequence 
of (3.7) and (3.8). Now suppose s(t;0, so) = s(t;0, so) where so,so G C+. 
Then (3.19) implies 

So - So = (1 - e"*) (<y5o(so) - <^o(so)) • 

Prom this and (3.31) we infer 

We observe that the derivative of this map is uniformly bounded away 
from zero. Indeed, (3.19) and (3.24) yield 

>e2*(l-K(so)|(l-e-*)) >e*. 

It follows by the inverse function theorem that fit is an open set, and by 
continuity the image of (C+ is fit. The inverse map from fit to C+ is analytic 
in fit, and on fit. 

Por ti > 0, the inverse of the map sq i— > si = s(ti; 0, sq) may be obtained 
by solving the characteristic equation in (3.16) backwards from time ti to 
to = 0, so we have so = s(0;ti,si). Now whenever ti > to > in general, 
we may follow any characteristic curve back from a point in fit^ at time ti 
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to a point in C_|_ at time and then forward to a point in at time to. 
This means that s{ti;tQ, sq) = s{ti;0, s{0;tQ, sq)). Part (iii) of the lemma 
now follows from the properties established in the case = 0. 

Next we prove part (i). For t > 0, Tf is the image of the map k i— > 
s(t;0,ik) = e^*{ik — ipQ(ik){l — e~*)), A; e M, and this is a function of k. 
We have s(t;0, 0) = 0, but Re s < for A; 7^ 0. This is so because Res and 
'Reipoiik) have opposite signs, and 

/•oo 

'R.e (faiik) = / {1 — coskx)nQ{x)dx > ^, A: / 0, 
Jo 

since no is continuous. Finally, we find that 

Im-^s(t;0,iA;) > e^\l - |uo(iA;)|(l - e"*)) > 

using (3.24). Hence Res is a function of Ims on F^. 

Now we establish part (ii). By (3.30) we have that as \sq\ 00 with 
So G C+, \^Pq{sq)/sq\ ^ 0, so s = soe2*(l + o(l)) by (3.19). Let si e C lie 
to the right of F^, and put /(so) = s(i;0, so) — si. It follows by applying 
the argument principle to large semicircles that the analytic function / has 
a single zero at some point so G C+. Indeed, arg/(re*^) ^ 9 as r ^ 00 
for — ^ < ^ < f , and as k goes from 00 to —00, f{ik) docs not cross the 
positive real axis, so argf{ik) changes from ^ to Thus, / maps a large 
semicircle to a curve that winds exactly once about 0. Hence si G Jlj. 

Finally, part (iv) follows by a change of variables, replacing t — to by t, 
and applying parts (i)-(iii). □ 



3.4 Proof of Theorem 3.1 

1. By the Fourier-Laplace inversion formula, it suffices to prove 



lim sup 

a;>0 



ikx 



v{t, ik) 



{ik)] dk 



0. 



(3.32) 



2. Let e G (0, |), and put R = ^e~^. We will prove (3.32) by estimating 
the integral for t > T separately in three regions: \k\ < R, R < \k\ < 
i?e^* and Re^^ < \k\, where R = Re^^^ and T depends only on e and 
the initial data vq. This is the same decomposition used in the proof of 
Theorem 2.1, and convergence in the region |A;| < i? will follow as before. 
However, estimates for |A:| > i? are more subtle, and use the analyticity and 
geometry of the characteristic map. 
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3. \k\ < R: Theorem 7.1 in [15] implies that xh{T,x)dx — xh*fi{x) dx 
in the sense of weak convergence of measures. Combined with (3.7) this also 
implies that the measures x'^h{T,x) dx converge weakly to x'^h^:^l{x) dx as 
t oo. But this implies v{t,ik) converges to uniformly on compact 

subsets of C+, and in particular on compact subsets of the imaginary axis [7, 
XV.3.2]. Thus, 

rR 

lim / \v{t, ik) — i{ik)\ dk = 0. (3.33) 

4- \k\ > R: It is sufficient to consider only k > R, since \v{t,ik)\ = 
\v{t, —ik)\. We will control v{t,ik) and separately: 

/■oo foo rco 

/ \v{t,ik) — v^^i{ik)\dk < / \v{t,ik)\dk+ / \v*^i{ik)\dk. 
Jr ' Jr Jr ' 

But \v^,iiik)\ < (2A;)-3/2 by (3.13). Thus, 



/"OO /"CXD 

/ \v^,i{ik)\dk < (2fc)-3/2 dfc = (2i?)- 
Jr ' Jr 



■1/2 



5. In the rest of the proof we estimate \v{t,ik)\dk. In order to aid 
the reader, we state the main estimates as two distinct lemmas. 

Lemma 3.4. Let e G (0, |). There exists T > depending on e and the 
initial data, and a universal constant C, such that if t>T then 



f 

Jr 



Re^(*-T) 

\v{t,ik)\dk <Cs. (3.34) 



Lemma 3.5. Let R > 0. There exists C depending on R and the initial 
data such that for allt>0 we have 

/•oo 

/ \v{t,ik)\dk <Ce-K (3.35) 

6. We now prove (3.32). We choose T as in Lemma 3.4, and then 
R = i?e~^^ in Lemma 3.5. Choose T* > T such that for t > 



R 

\v{t, ik) - v*^i{ik)\ dk < e, (7e~* < Ce"'^* < e. 

R 
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Thus, for t > T* we have 

/ \v{t,ik) — v*^i{ik)\dk < / \v{t,ik) — v*^i{ik)\dk 
Jr ' J-R 

(roo f'Re^^ poo 

/ \v^i{ik)\dk+ / \v{t,ik)\dk+ / \v{t,ik)\dk 
Jr ' Jr JRe'^t 

<£ + 2{e + C£ + e). 
Since £ G (0) |) niay be chosen arbitrarily small, this completes the proof. 

3.5 Proof of Lemma 3.4 

In this subsection we will always suppose s G C+. In a manner similar to step 
6 of the proof of Theorem 2.1, the idea is to get estimates on the semicircle 

Cr '■= {s G C+ I |s| = R} valid for large time, and propagate these estimates 
outwards along characteristics. We first use (3.10) and (3.13) to obtain the 
following estimates for s G C+: 

|^*,i(s)| < |2s|V2, |u*,i(s)| < |2s|-V2, |^;^^^(s)| < |2s|-3/2. (3.36) 

Next, we use the uniform convergence on compact sets and (3.36) to see 
that there exists Tq (depending on e and the initial data) such that for all 
So G Cr and to > Tq we have 

\ip{to,so)/so\ < 2(2i?)-V2 = 2£ < 1/4, (3.37) 
|n(to,so)| < i2R)-^/^=e, (3.38) 
Hto,so)\ < (2i?)-3/2 = £3_ (339) 

We first extend (3.37) to a larger domain in s. 
Claim 1: There exists Ti > Tq such that 

< 1/3, t>Ti, s G C+, |s| > R. (3.40) 

Proof of claim 1: Observe that by using (3.27) and (3.30) in (3.21), we have 

a := sup{|</?(To, s)/s\ \ s G C+, |s| > ii!} < 1. 

Fix ti > Tq, Si G C+ with |si| > R. Either the characteristic curve s{t; ti, si) 
that passes through si at time ti intersects Cr at some time to G [To,ii], 



^{t, s) 
s 
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or not. If so, then si = s{ti]to,so) for some sq G Cr, and (3.21) and (3.37) 
directly yield 



(p{tl,Sl] 



Sl 



< 



1/4 



1 - 1/4 



1 

3' 



If not, then |s(t;ii,si)| > R for all t G [ro,ii], by continuity and the fact 
that s{t;ti,si) G C+ for all t G [0,ti] by part (iv) of Lemma 3.3. Then 
taking to = Tq, sq = ^(ro;*!, si) in (3.21) yields 



^g-(tl-To) I 
< < -, 

- 1-a -3' 



provided ti > Ti with Ti sufficiently large. This proves the claim. 
Claim 2: Let T = Ti + iln2. Suppose ii > T and < |si| < Re^^*^-'^\ 
Then the characteristic curve s{t;ti,si) that passes through si at time ti 
intersects Cr at some time to S [Ti,ti]. 

Proof of claim 2:. Suppose the claim were false. Then the continuity 
of |s(t;ti,si)| and part (iv) of Lemma 3.3 imply R < \s{to;ti, si)\ for all 
to G [Tijti]. But now, by (3.20) with sq = s(to;ti,si) we have 



so 1- 



V^ito, So) 



So 



-2(ti-to). 



Sl 1- 



(p{ti,si) 



Sl 



(3.41) 



We take to = ^i aiid apply (3.40) and the hypothesis < Re^^^^ = 
^i?e2(*i-^i) to deduce 

i?<|.o|<|si|e-2(*-^^)^^<i?, 

a contradiction. This proves the claim. 

We now apply these claims to propagate the decay estimate (3.39). Prom 
claim 2, for any t = ti>T,R<k< i?e^'^*~^\ with si = ik we obtain 
*o e [^1,*] and Sq G Cr and substitute (3.20), (3.39) and (3.40) in the decay 
estimate (3.26) to obtain 



\v{t,ik)\ < 



\v{to,so)\ 



(1- |n(to,so)|)3 



So - (p{to,so) 



ik — ip{t, ik) 



3/2 



< (l-£)-3 |i;(to,so)| 



2so 



3/2 



< (l-£)-3(2i?)-3/2 (^^y^ = {l-er'k-^/\ 
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Therefore, 

/ \v(t,ik)\dk<(l-e)-^ k'^''^ dk = < 

Jr Jr (1 - 

with C = 2(8/7)^2^2. This completes the proof of Lemma 3.4. 



3.6 Proof of Lemma 3.5 

We consider the initial time to = and the following special case of (3.19): 
s = s{t; 0, so) = e^* [sq - m{so){l - e"*)] . (3.43) 
For any t>0, the map sq h- >■ s{t; 0, Sq) is analytic for Re(so) > 0, and 

ds 

— = e^' (1 - txo(so)(l - e-*)) , uo(so) = u(0, so). (3.44) 

Recall that denotes the image of C+ under sq ^ s{t] 0, sq), and Tf denotes 
the image of the imaginary axis; we let T^t denote its preimage. As was 
observed in Lemma 3.3, Tf is a graph over the imaginary axis, contained in 
the left half plane. 

We will use the analyticity of v{t, s) in and contour deformation. For 
large finite R2 < oo consider the domain ABCD shown in Figure 3.6. The 
path AB is chosen so that A'B' is a straight line. CD is parallel to the real 
axis, and lies in 0^ since Ft is a graph over the imaginary axis. Then by 
Cauchy's theorem, 

rR2 r 

/ e'^''v{t,ik)dk= / e'^''v{t,ik)dk 

jRe^t JbC 

= / e^''v{t,s)ds + / e^''v{t,s)ds + / e^^'vit, s) ds. 

JDA JAB JCD 

Let a denote Re s. Since cr < in J7t for s G CD we see that the last integral 
is estimated by 

[ e'-v{t,s)ds < sup \v{t,s)\ r e-da=^^^P26CZ^^!Ml. 

JcD s&CD J-00 X 

By the decay estimate (3.26) we have 



\voiso)\e 3* 

seCD ' ' ' " ~ sieCD (1 - |uo(so)|) 



sup \v{t,s)\< sup — ; — - — So = s{0;t,si). 
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Figure 3.1: The s-planc is on the left, the .sp-plane on the right, i^t is the 
region to the right of Tt. A = s{t; 0, iR), B = iRe^\ C = ii?2, Im(D) = R2, 
A' = ik, B' = s{0;t,iRe'^^) 



It follows from (3.30) and the fact that |so| = l^ile ^*(1 + o(l)) — 00 as 
R2 ^ 00 that supg^^QD |t;o(so)| 0- We thus let i?2 — 00 to conclude that 



poo f r 

/ e'^''v{t,ik)dk= / e'''v{t,s)ds + / e'''v{t,s)ds. 



(3.45) 



where T-t^A denotes the path from 00 to A on Fj. Notice that (3.45) holds 
independent of x. 

The virtue of deforming the contour is that the integrals may now be 

estimated by changing variables from s to sq. We use the solution formula 
(3.23) together with the change of variables s = s{t;0,ik) and (3.44) to 
obtain 



[ e^^v{t, s)ds = ie * j e 

JVtA JR 



s{t;0,ik)x _ 



voiik) 



(1 - uo{ik)il - e-^)y 



■dk. 



Since Res(i;0,zfc) < and sup,i^,^pg^\uQ{ik)\ < 1, this yields the estimate 



e^^v{t, s) ds 



< Cie ^\\vo\\li. 



(3.46) 
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Similarly, we have by (3.23) and (3.44) 

/ e'''v{t,s)ds = e"* / e^(*'°'^o)^- 
Jab Ja'b' (1 



vo{so) 



■ dst 



{l-uo{so){l-e-*)y 

< e-^\A'B'\ sup |l-?xo(so)(l-e-*)^^ 
soeA'B' 

The point A' = iR is independent of t. It also follows from (3.43) that B' = 
s{0;t,iRe'^^) converges to the point sq G C+ that solves iR = sq — ipo{so)- 
Thus, we have the exponential decay estimate 



/ ■ 

Jab 



e''''v{t, s) ds 



< ^26-*. 



(3.47) 



The constants in (3.46) and (3.47) depend only on R and the initial data 
uq. To be explicit, we set C = Ci\\vo\\li + €2- This completes the proof. 



3.7 The discrete Smoluchowski equations 

We now use the proof of Theorem 3.1 to obtain a uniform convergence 
theorem for the discrete Smoluchowski equations with additive kernel. The 
proof is simpler and we do not need the contour deformation argument. 

Let Vt = 'Yl^i''^li^)^hl{x) denote a measure-valued solution to (1.1). We 
first adapt the rescaling (3.4) and (3.5) to similarity variables. Let 

l = lhe-'^\ ni{t) = h-^e^*ni{t). (3.48) 

Then the discrete Fourier inversion formula analogous to (2.32) is 



PMt) = 5^ / 



e'%{t, ik) dk. 



(3.49) 



Theorem 3.6. Let i'q > be a lattice measure with span h such that 
xi'o{dx) = x'^i'o{dx) = 1. Then with the scaling (3.4-8) we have 



lim sup I 



n 



ni{t) - n* i(Z) 



= 0. 



Proof. By (3.49) and the continuous Fourier inversion formulas it suffices to 
show that 



lim sup 



e'^''''v{t,ik)dk- I e'""'v^,i{ik)dk 



I 

Jr 



-ilhk„ 



(3.50) 
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Let £ G (0, |) and choose R = ^e"^- The integrals over [—R,R] and R < 
\k\ < i?,e^* with R = e~'^^ arc controlled as in the proof of Theorem 3.1, and 
it only remains to control the integral of \v{t,ik)\ over i?e^* < k < 7re^*//i. 
This is considerably simpler than in the previous proof. We use the solution 
formula (3.23) and change variables via ik = s{t; 0, sq), using (3.44) to obtain 



e''^''v{t, ik) dk = ie 



-t 



r_t(R,7r//i) (1 - uo{so){l - e-*))' 



dsQ. 



Here r_t(i?, vr/Zi) denotes the segment along the curve T t from s(0; t, iRe^*) 
to s{Q;t,me^^ /h). The formula (3.19) shows that r_t(^,7r//i) converges to 
a compact curve defined implicitly by ik = sq — ipo{so),R < k < ir/h. 
Thus, for t > r we have 



/ 

7r_t(R,7r//i) (1 



■ dso 



uo(so)(l - e-*)) 
Thus, this term is less than e for all t large enough. 



<C{T,R,uo,vo)e 



-t 



□ 



4 Self-similar gelation for the multiplicative kernel 

For K = xy, McLeod solved the coagulation equation explicitly for monodis- 
perse initial data, and showed that a mass-conserving solution failed to exist 
for t > 1. The second moment satisfies m2{t) = (1 — The divergence of 
the second moment indicates that breakdown is associated with an explosive 
flux of mass toward large clusters. A rescaled limit of McLeod's solution is 
the following self-similar solution for K = xy [1]: 

n(t,x) = ^a;-5/2e-(l-*)'^/^ x > 0, t < 1. (4.1) 

Evidently this solution has infinite mass (first moment) . This should not be 
thought unnatural, however, since it was shown in [15] that equation (1.1) 
has a unique weak solution for any initial distribution with finite second 
moment. 

The problem of solving Smoluchowski's equation with multiplicative ker- 
nel can be reduced to that for the additive kernel by a change of variables [4]. 
Let us briefly review this. In unsealed variables we define 

POO 

^it,z)= {l-e-''')xn{t,x)dx. (4.2) 
Jo 
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Then * solves the inviscid Burgers equation: 

- ^5^* = 0, (4.3) 

with initial data 

^o(^)=/ (1 -e-^'^)xno(x)dx. (4.4) 

JO 

The gelation time for initial data with finite second moment is Tgei = 
{j^ x'^VQ{dx))~^ and this is exactly the time for the first intersection of 
characteristics [15]. We presume that the initial data is scaled to ensure 

roo poo 

/ x'^nQ{x)dx= / x^nQ{x) dx = 1. (4.5) 
Jo Jo 

Then the gelation time is Tgd = 1. The connection between the additive 
and multiplicative kernels is that ^' solves (4.3) with initial data ^'o, if and 
only if ^{t,z) is a solution to (3.3) with the same initial data, where 

^{t, z) = p7$(r, z), with r = log(l - t)'^ . (4.6) 

For solutions n°'"^(t, x) and n^^ij, x) to Smoluchowski's equation with mul- 
tiplicative and additive kernels respectively, this means that 

XTr'\t, x) = {\- t)-^n^^^(r, x) (4.7) 

for all t G (0, 1), if and only if the same holds at t = 0. We thus obtain a 
scaling limit as t Tgd directly from Theorem 3.1. The similarity variables 
for the multiplicative kernel are 

n9 ^, n(t, x(l — n(t,x) , 

X=il-tfx, hit,x)= (i_^^5 -(ff^' (4-8) 

and the self-similar profile is 

n.,2(x) = -T^e-'/'. (4.9) 

Theorem 4.1. Suppose no{x) > 0, J^x^no{x)dx = x^no{x) dx = 1. 
Suppose also that the Fourier transform ofx^riQ is integrable. Then in terms 
of the rescaling (4-8) we have 

lim sup x^\h{t,x) — h^^2{x)\ = 0, (4-10) 
where n*,2(^^) is the self-similar density in (4-9). 
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Notice that (4.8) is not a mass-preserving rescaling; indeed, the rescaled 
mass diverges: 



> oo. 
1-t 



I xn{t,x)dx = - / xn{t,x)dx = - — 

Jo ^ — t Jo 1 ~ 

Instead, (4.8) preserves the second moment: 

f'co poo 

/ Sp'h{t,x)dx = il-t) j x'^n{t,x)dx = 1, te[0,l). 
Jo Jo 

The explanation is that the scaling in (4.8) is designed to capture the behav- 
ior of the distribution of large clusters as t approaches Tg^i — the average 
cluster size is {l—t)~^. Correspondingly, the mass of the self-similar solution 
is infinite. 

Theorem 3.6 may be similarly adapted to K = xy. In the discrete 
case, the correspondence (4.7) between solutions of Smoluchowski's equa- 
tions with multiplicative and additive kernels becomes 

hlnr\t) = (1 - ty'nf ''{logil - t)-') (4.11) 

We introduce similarity variables via 

l = lh{l-tf, ni{t) = h~\l-t)-^ni{t). (4.12) 

Then directly from Theorem 3.6 we obtain the following. 

Theorem 4.2. Let i^q > be a lattiee measure with span h such that 
x'^i'o{dx) = x^i'o{dx) = 1. Then with the rescaling (4-12) we have 



lim sup I 



ri 



Mt) - fi*,2{i) =0. (4.13) 
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